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Abstract 

I discuss the spherically symmetric but inhomogeneous Lemaitre-Tolman- 
Bondi (LTB) metric, which provides an exact toy model for an inhomogeneous 
universe. Since we observe light rays from the past light cone, not the expansion 
of the universe, spatial variation in matter density and Hubble rate can have the 
same effect on redshift as acceleration in a perfectly homogeneous universe. As 
a consequence, a simple spatial variation in the Hubble rate can account for the 
distant supernova data in a dust universe without any dark energy. I also review 
various attempts towards a semirealistic description of the universe based on the 
LTB model. 
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1 Introduction 



The simplest homogeneous and isotropic cosmological models, based on the Friedmann- 
Robertson-Walker (FRW) metric, have proved to be remarkably successful ever since 
Edwin Hubble in 1929 rather cautiously suggested that the apparent linear correlation 
between the observed redshifts and distances of 24 galaxies could hint towards the 
possibility "that the velocity-distance relation may represent the de Sitter effect, and 
hence that numerical data may be introduced into discussions of the general curvature 
of space" flU. Indeed, numerical data now guides the development of cosmology, 
which has become a precision science, albeit mostly within the framework of a 
perfectly homogenous background metric. 

The FRW universe is characterized by two functions, the Hubble rate H and 
the density parameter fi, or the average expansion rate and the average density of 
mass energy, respectively, which depend on time but are independent of the spatial 
location. However, one should keep in mind that their values cannot be extracted 
directly from the observations but must be deduced from the properties of light coming 
from the past light cone. In the context of the FRW model this is almost trivial, since 
the redshift z and scale factor a(t) are everywhere related by z = a(t )/a(t e ) — 1, 
where the subscripts refer respectively to the observation and the emission of light. 
This theoretical simplicity should however not cloud the fact that all cosmological 
parameter determination requires an element of interpretation of the data. Of course, 
the FRW interpretation of the properties of the past light cone has served cosmology 
well, giving a good fit to observations and, until the late 90's, implying a matter 
dominated universe with £1 m £l M . 

The situation changed dramatically with the WMAP [2] and distant supernova data 
0. Considering the recent data from supernovae EMU, galaxy distributions [6] and 
anisotropics of the cosmic microwave background [7], the simplest FRW model would 
now lead to a highly contradictory picture of the universe, with the following best fit 
values for the average matter density: 

• Cosmic microwave background: £l M ~ 1 

• Galaxy surveys: Hm ~ 0.3 

• Type la supernovae: Qm ~ 

As is well known, the glaring discrepancies between the different data sets have 
conventionally been remedied by introducing the cosmological constant A or vacuum 
energy £^a to the Einstein equations. This gives rise to an accelerated expansion 
of the universe. As a consequence, the apparent dimming of the luminosity of 
distant supernovae finds, in the context of perfectly homogeneous universe, a natural 
explanation 1 . 

However, although the cosmological concordance ACDM-model [H fits the ob- 
servations well, there is no theoretical understanding of the origin of the cosmological 

'Even if the primordial perturbation is not scale free, the combination of the CMB fluctuations 
and the shape of the correlation function up to ~ 100/i _1 Mpc, seems to require dark energy for a 
homogeneous FRW model [ 10 1. 
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constant or its magnitude. For particles physicists, who have spent a long time 
trying to prove that the cosmological constant must be zero, the tremendously small 
cosmological constant which just now happens to start to dominate the energy budget 
of the universe, is a theoretical nightmare. There exist a large number of different 
dark energy models (see e.g. (Hill) that attempt to provide a dynamical explanation 
for the cosmological constant, but none of them are compelling from particle physics 
point of view; moreover, very often they require fine-tuning. Modifications of the 
general theory of relativity on cosmological scales appear to suffer from analogous 
problems. For instance, f(R) gravity theories |[TT1 in the metric formalism are plagued 
by instabilities Ifl2l while in the Palatini approach the cosmological constant seems to 
be essentially the only consistent modification that fits all the cosmological data [TT3ll . 

Facing such difficulties, one might be tempted to consider relinquishing the FRW 
assumption of the perfect homogeneity of the universe. After all, inhomogeneities 
are abundant in the universe: there are not only clusters of galaxies but also large 
voids. Because general relativity is a non-linear theory, even relatively small local 
inhomogeneities with a sufficiently large density contrast could in principle give rise 
to cosmological evolution that is not accessed by the usual cosmological perturbation 
theory in an FRW background. In fact, the potentially interesting consequences of 
the inhomogeneities were recognized already at the time when the homogeneous and 
isotropic models of the universe were first studied, but their impact on the global 
dynamics of the universe is still largely unknown (see e.g. O). Then the question 
arises: could the acceleration of the universe be just a trick of light, a misinterpretation 
that arises due to the oversimplification of the real, inhomogeneous universe inherent 
in the FRW model? Light, while traveling though inhomogeneities, does not see 
the average Hubble expansion but rather feels its variations, which could sum up to 
an important correction 2 . This effect is particularly important for the case of large 
scale inhomogeneities which will be the focus of the present paper. If the local 
Hubble expansion rate were to vary smoothly at scales of the order of, say, thousand 
megaparsecs, that would very much change our interpretation of the distant supernova 
redshifts. In such an inhomogeneous universe we could also just happen to be located 
in a special position. For instance, fate could have relegated us to an underdense region 
with a larger than average local Hubble parameter so that the discrepancy between 
nearby and distant supernovae luminosities could be resolved without dark energy. 

Local inhomogeneities have recently been invoked as the culprit for the apparent 
acceleration of the expansion of the universe 3 , in particular by virtue of their so- 
called backreaction on the metric (for a discussion on the issues involved and a 
comprehensive list of references, see lfT71[T8lO . One constructs an effective description 
of the universe by averaging out the inhomogeneities to obtain averaged, effective 
Einstein equations which, in addition to the terms found in the usual homogeneous 
case, include new terms that represent the effect of the inhomogeneities IU91I201I2TTI . 

However, since we can only observe the redshift and energy flux of light arriving 

2 For a recent calculation of the small scale inhomogeneity-induced correction to the cosmological 
constant that one would infer from an analysis of the luminosities and redshifts of Type la supernovae, 
assuming a homogeneous universe, see 03]. 

3 Inhomogeneities as an alternative to dark energy were first discussed in 1 16 1. 
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from a given source, not the expansion rate or the matter density of the universe 
nor their averages, one may wonder what are the actual observables related to the 
averaged equations. To wit, since we do not observe the average expansion of 
the universe directly, its average acceleration is also an indirect conclusion, arising 
from the fact that in the perfectly homogeneous cosmological models dark energy is 
required for a good fit. Consequently, there is no a priori reason to assume that an 
accelerated expansion is necessarily required to fit the data if one assumes a general 
inhomogeneous model of the universe. One may also add that the averaging procedure 
as such is not without problems: in general it is not correct to integrate out constrained 
degrees of freedom as if they were independent, and in cosmology the fact that we 
can make observations only along our past light cone makes the observable universe 
a constrained system. Hence it would be desirable to study the effects of the the 
inhomogeneities on the directly observable light in an exact cosmological model. 
Unfortunately, in the presence of generic inhomogeneities this would be practically 
an impossible task. Instead, one must resort to toy models, the simplest of which is 
the spherically symmetric but inhomogeneous Lemaitre-Tolman-Bondi (LTB) model 
J22l|23ll2l. ' 

The great virtue of the LTB model is that it is exact. Because of its high degree of 
symmetry, it may not be realistic as such, but the LTB model is nevertheless interesting 
at least on two counts. First, it serves as a simple testing ground for the effects of 
inhomogeneities when fitting the cosmological data without dark energy. Second, since 
the fits can be performed unambiguously, the nature of the effective accelaration in the 
models where the spatial degrees of freedom have been averaged out, can be made 
transparent by comparing the averaged and "exact" models. 

Of course, one can also take the LTB model more seriously. For instance, one may 
use the LTB metric to describe a local underdense bubble in FRW universe, for which 
there is some evidence both from supernova [25 J and galaxy data E6l . First attempts 
along these directions ETll assumed an underdense region separated from the outside 
homogeneous FRW universe by a singular mass shell, followed by investigations of 
more realistic models with a continuous transition between the inner underdensity and 
the outer homogeneous universe (see e.g. Il28l |29l ). More complicated situations, 
including off-centered observers, can also be addressed, as will be discussed in Sect. 

m 

2 The Lemaitre-Tolman-Bondi metric 

Let us consider a spherically symmetric dust universe with radial inhomogeneities as 
seen from our location at the center. Choosing spatial coordinates to comove (dx % / dt = 
0) with the matter, the spatial origin (x l = 0) as the symmetry center, and the time 
coordinate (a; = t) to measure the proper time of the comoving fluid, the line element 
takes the general form 1221 1251 1251 

ds 2 = -dt 2 + X 2 (r, t)dr 2 + A 2 (r, t) (d6 2 + sin 2 9d<p 2 ) , (2.1) 
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where the functions A(r,t) and X(r,t) have both temporal and spatial dependence. 
The homogeneous FRW-metric is a special case and is obtained by letting 

X(r,t)^-^=, A(r,t) -> a(t)r. (2.2) 

The energy momentum tensor is given by 

T^ = -p M (r,t)8%5°„-p A 8%, (2.3) 

where p M (r,t) is the matter density, = 5q represent the components of the 4- 
velocity-field of the fluid, and we have kept the vacuum energy p A for generality. Note 
that although the fluid is staying at fixed spatial coordinates, it can physically move in 
the radial direction. Plugging Eq. (12.11) into the Einstein equation, G^ u = 8irGT%, 
one finds the set of equations 



87tG(p m + pa), (2.4) 
(2.5) 



and 




8nGp A , (2.6) 



A" A AX A'X' X 

+ T + T^ + -r^T + ^7 = 8nGp A . (2.7) 



AX 2 A AX AX 3 X 



These contain only three independent differential equations, and we may solve X and 
X from Eq. (12.51) and A' 2 and A" from Eq. (12.61) . Then one can substitute these into 
Eq. (12.71) and find that it yields an identity. Thus only two of equations (I2.5I) - (I2.7I) are 
independent. One can easily solve Eq. (12.51) to obtain 



X(r,t) = C(r)A'(r,t) , (2.8) 
where the function C(r) depends only on the coordinate r. By redefining C(r) = 



1 / a/1 — k{r), where k{r) < 1, we can thus write the LTB metric Eq. (12.11) in its usual 
form: 

ds 2 = -dt 2 + ^'^P' dr 2 + A 2 (r, t) (d9 2 + sin 2 9d^ 2 ) , (2.9) 
1 — k(r) x ' 

where k(r) is a function associated with the curvature of t = const, hypersurfaces. 
The FRW metric is the limit A(r, t) — > a(t)r and k{r) — > kr 2 . 
The two independent equations are given by 

A 2 + k(r) 2AA' + k'(r) 

+ AA , = 8vrG(p A/ + PA), (2.10) 

A 2 + 2AA + k(r) = 8nGp A A 2 . (2.1 1) 
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The first integral of Eq. (12.111) is 



A 2 



F(r) 8ttG k(r) 
-JT + — Pa - 



(2.12) 



where F(r) is a non-negative function that, like k(r), is fixed by the boundary 
condition. Substituting Eq. (|2.12l) into Eq. (|2.10l) yields 



F' 
A' A 2 



MGp 



M 



(2.13) 



By combining Eqs. (12.101) and (12.1 II) we can construct the generalized acceleration 
equation 

2 A I A' AixG , 

ZA + 3A> = ~— (P " ~ ^ (2 ' 14) 
which implies that the total acceleration, represented by the left hand side, is negative 
everywhere unless the vacuum energy is large enough: p A > pm/2. However, it does 
not exclude the possibility of having radial acceleration (A'(r, t) > 0), even in the pure 
dust universe, if the angular scale factor A(r, t) is decelerating fast enough, and vice 
versa. This serves to demonstrate how the very notion of the acceleration becomes 
ambiguous in the presence of the inhomogeneities [|30l . 

The boundary condition functions F(r) and k{r) are specified by the exact 
physical nature of the inhomogeneities. Their relation to the FRW model parameters 
can be recognized by comparing Eq. (12.121) with the Einstein equation for the 
homogeneous FRW-model 



H\t) = 



a(t) 
a(t) 



8nG / . k 

— (Pm + Pa) - -J 



Hi 



n 



M 



a J 



(2.15) 
(2.16) 



where a = a(t ) and H = H(t ). Thus, a comparison between Eqs. (12.121) and 
(12.151) motivates one to define the local Hubble rate as 



H(r,t) 



A{r,t) 



(2.17) 



A(r,t) • 

The local matter density can be defined through 

F(r) = H 2 Q (r)n M (r)A 3 (r) , (2.18) 

with 

k(r) = H 2 (r)(n M (r) + Q A (r) - l)A 2 (r) , (2.19) 

where we have defined the boundary values at t through A (r) = A(r, t ), H (r) = 
H(r, t ), and Oa(?") = 8nGpA/SHl(r). With these definitions, the position-dependent 
Hubble rate, Eq. (12.121) . takes a physically transparent form OTA : 



H 2 (r,t)=H 2 Q (r) 



A 



A 



(2.20) 
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where fl c (r) = 1 — £Ia(t) — Om(f). 

The difference between the conventional Friedmann equation (12.151) and its LTB 
generalization, Eq. (12.201) . is that all the quantities in the LTB case depend on the r- 
coordinate. Thus in the presence of inhomogeneities, the values of the Hubble rate 
and the matter density can vary at every spatial point so that the inhomogeneous 
dust models are defined by two functions of the spatial coordinates: H (x z ) and 
VLm{x 1 ). As a consequence, the inhomogeneities are of two physically different kinds: 
inhomogeneities in the matter distribution, and inhomogeneities in the expansion 
rate. Although their dynamics are coupled via the Einstein equation, as boundary 
conditions they are independent. The universe could have an inhomogeneous big bang, 
where the universe came into being at different times at different points, and/or an 
inhomogeneous matter density. This opens up the possibility for an inhomogeneous 
universe that has a homogeneous present-day VLm\ a model of this kind could 
potentially fit the supernova data as well as the galaxy surveys without invoking dark 
energy. However, if VL M {r) = const., the physical matter distribution p M itself has 
a spatial dependence provided H (r) ^ const.. It can be made constant by choosing 
VL M {r)Hl{r) = const. 

The spatial dependence holds true even for the gauge freedom of the scale function. 
In the FRW case the present value of the scale factor a(t ) can be chosen to be any 
positive number. Similarly, the corresponding present-day scale function A(r,t ) of 
the LTB model can be chosen to be any smooth and invertible positive function. In 
what follows we will choose the conventional gauge 



Integrating Eq. (12.201) then gives the relation between the scale factor A(r, t) and the 
coordinates r and t, which can also be used to find the age of the LTB universe. One 
finds 



For any space-time point with coordinates (t, r, 9, ip), Eq. (12.221) determines the 
function A(r, t) and all its derivatives. Thus the metric Eq. (12.91 ) is specified, and given 
the inhomogeneities, all the observable quantities can be computed. Eq. (|2.22l) can be 
integrated in terms of elementary functions when £l\(r) = or f^r) + fW(r) = 1; 
as an example, in the latter case one finds 



A(r,t ) 



= r . 



(2.21) 




(2.22) 



(t-t )H { 



o — 



3y/l-n M (r) 



2 




arsinh-y/ u(r) , (2.23) 



where 




(2.24) 



In this particular case A(r, t) can be found explicitly as 




(2.25) 
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where r = ^/QnGp\{t — t ). 



3 Inhomogeneities and luminosity distance 

To compare the inhomogeneous LTB model e.g. with the supernova observations, 
we need an equation that relates the redshift and energy flux of light with the exact 
nature of the inhomogeneities. For this, one must study propagation of light in the 
LTB universe 4 . Let us here derive the appropriate equations for notational clarity; a 
more general derivation for an off-center observer can be found in ll33l . 

From the symmetry of the situation, it is clear that light can travel radially, that is, 
there exist geodesies with d6 = dip = 0. Moreover, since light always travels along 
null geodesies, we have ds 2 = 0. Inserting these conditions into the equation for the 
line element, Eq. (|2.9I) . we obtain the constraint equation for light rays 

dt_ = dr A'(r,t) 
du du y/l - k(r) ' 

where u is a curve parameter, and the minus sign indicates that we are studying radially 
incoming light rays. 

Consider two light rays with solutions to Eq. (13.11) given by t\ = t(u) and t 2 = 
t(u) + X(u). Inserting these to Eq. (|3.1I) we obtain 

d dt(u) dr A'(r,t) 

-i-h = — ^ = -- . v ' ' (3.2) 

du du du y/i — k(r) 

d t _ dt(u) + d\(u) _ dr A'(r,t) + dX(u) 
du du du du «/l — k(r) du 

d t _ dr A'(r,t(u) + X(u)) _ dr A'(r,t) + A'(r,t)X(u) 
du du y/l - k(r) du ^/l - k(r) 

where Taylor expansion has been used in the last step and only terms linear in X(u) 
have been kept. Combining the right hand sides of Eqs. (13.31) and (13.41) gives the 
equality 

dX(u) dr A'(r,t)X(u) 



du du yj\ - k(r) 

Differentiating the definition of the redshift, z = (A(0) — X(u))/X(u), we obtain 

dz _ dX(u) A(0) _ dr (1 + z)A'(r,t) 
du du X 2 (u) du y/l — k(r) 



(3.5) 



(3.6) 



Luminosity distance in a perturbed FRW universe has been considered in [32] 
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where in the last step we have used Eq. (13.51) and the definition of the redshift. Finally, 
we can combine Eqs. (12.191) . (13.11) and (13.61) to obtain the pair of differential equations 

dt 

dz 

dr 
dz 

which determine the relations between the coordinates and the observable redshift, i.e. 

t(z) and r{z). 

Now that we have related the redshift to the inhomogeneities, we still need the 
relation between the redshift and the energy flux F, or the luminosity-distance, defined 
as d L = ^L/AtcF, where L is the total power radiated by the source. This is given by 

m 

d L {z) = {l + z) 2 A{r{z),t{z)). (3.9) 
Likewise, the angular distance diameter is given by 

d A (z) =A(r(z),t(z)) . (3.10) 

As the ^-dependence of t and r are determined by Eqs. (|3.7I) and (13.81) and the scale 
function A(r, t) by Eq. (12.221) . using Eq. (13.91) one can calculate g?l for a given z. 
All of these relations have a manifest dependence on the inhomogeneities (i.e. on the 
functions H (r) and Q M (r)). What remains is a comparison of Eq. (13.91) with the 
observed d L (z). 

Because the boundary functions of the LTB model are arbitrary, it comes as no 
surprise that any isotropic set of observations can be explained by the appropriate 
inhomogeneities of the LTB model 11351 . That the supernova data could be interpreted 
in terms of an inhomogeneous LTB model with no cosmological constant was first 
suggested by Celerier 061 . who pointed out that the LTB model is degenerate with 
respect to any magnitude-redshift relation so that the accelerated expansion could be 
modelled by a very large number of inhomogeneity profiles. In this sense the LTB 
model is not predictive. The intriguing aspect here is rather the matter of principle 
which the LTB model can be used to demonstrate: that the supernova data does not 
necessarily imply accelerating expansion and hence the existence of dark energy is 
not an unavoidable consequence of the data but rather depends on the framework the 
data is interpreted in. Moreover, the inhomogeneities need not contradict the observed 
homogeneity in galaxy surveys 0, as is often claimed (see e.g. 071 ). since the model 
admits solutions with constant tt M but with a position-dependent H. 

To demonstrate this, let us consider the gold sample of 157 supernovae of Riess 
et. al. [4] and disregard LSS and CMB data for the moment. In the FRW model the 
parameters that best describe our universe are found by maximizing the likelihood 
function exp(— x 2 (H , VLm, ^a)) constructed from the observations. However, to find 
the boundary conditions of the LTB universe that best describe our universe, we should 
in principle maximize the likelihood functional exp(— x 2 (H (r), n M (r))). In practice, 
this is impossible. One can only consider some physically motivated types for the 
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-A'(r,t) 
(l + z)A'(r,t) 

a/i + ^o 2 W(i 



n M (r)-n A (r))Al(r) 



(l + z)A'(r,t) 



(3.7) 
(3.8) 



functions H (r) and QmW that contain free parameters; these are then fitted to 
the supernova observations by maximizing the leftover likelihood function. In the 
literature there exist several fits to the supernova data employing a simple LTB model 
with different authors having chosen different density profiles (and, unfortunately, 
often a different notation) [|38ll28l[39ll40l[3Tll4ll. 

Since the expansion rate of the FRW universe has to accelerate in order to fit 
the supernova data, the second time derivative of the FRW scale function should be 
positive. In contrast, in the LTB universe the observations are affected by the variation 
of all the dynamical quantities along the past light cone, not just the time variation. 
Indeed, the directional derivative along the past light cone is given by 

d d dt d d A'(r, t) d d d 

— = 1 = , — ~ , (3-11) 

dt dt dr dr dt */l — k(r) 9r dt dr ' 

where the approximation in the last step is more accurate for small values of r, but is 
qualitatively correct even for larger r. 

The main message of Eq. (13.1 II) is that from the observational point of view, 
the negative r-derivative roughly corresponds to the positive time derivative. This is 
natural since by looking at a source, we simultaneously look into the past (i.e. along 
the negative t-axis) and into a spatial distance (i.e. along the positive r-axis). Hence, 
to mimic the acceleration, i.e. for the expansion rate to look as if it were to increase 
towards us along the past light cone, the expansion H (r) must decrease as r grows: 
hence we should look for an LTB model with H' (r) < 0. Thus, keeping in mind the 
homogeneity of galaxy distributions, we could choose a simple four parameter LTB 
model like OTTl 

H (r) = H + AHe~ r/r °, 

^m(^) = = constant , (3.12) 

where H, AH, Q an d r are f ree parameters determined by the supernova observa- 
tions. The best fit values are found to be OTTi 

H + AH = 66.8 km/s/Mpc, AH = 10.5 km/s/Mpc, r = 500 Mpc, Q = 0.45 . 

(3.13) 

The goodness of the fit is x 2 — 172.6 (x 2 /157 = 1.10). The confidence level 
contours with fixed values of f2o and H are shown in Fig. 1. For comparison with 
the homogeneous case, the best fit nonflat ACDM has VLm = 0.5, il\ = 1.0 with 
X 2 = 175 (x 2 /157 = 1.11). What is perhaps surprising is the fact that the supernova 
fit is not only in qualitative agreement with the observed homogeneity in galaxy 
surveys but also automatically yields a value for the present-day matter density that 
is consistent with the observations. The smallness (~ 15%) of the spatial variation in 
the Hubble parameter is also somewhat surprising, considering that it is of the same 
order as the uncertainty of the model-independent 5 determination of the local Hubble 
rate by the Hubble space telescope [|42l . The variation of the Hubble parameter found 

5 Note that the smaller uncertainties found in the CMB data analysis cannot be used here as those fits 
assume that the entire universe is perturbatively close to the homogeneous FRW model. 
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Figure 1: Confidence level contours in the LTB model with f2A/( r ) = constant = 0.45 
and H (r) = 56.3 km/s/Mpc + AHe~ r/r °. From EQ. 
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by Alnes, Amarzguioui and Gr0n [|28l . who used a different model Ansatz, has also 
similar magnitude, but in contrast their model contains a large (~ 400%) variation in 
the matter density at scales larger than the current range of galaxy surveys. 

One can also fit data with Eq. (13.121) together with a cosmological constant. Taking 
^M( r ) + ^A( r ) = 1 one finds no improvement PTTl . If instead of Q, M = const, we 
assume a strictly uniform present-day matter distribution with pu(r, t ) = constant, 
which implies Hq (r)ttM(r) = constant, we may choose the parametrization 

H (r) = H + AHe~ r/r °, 

n M (r) = n (H + AH) 2 /(H + AHe' r ' r °) 2 . (3.14) 

The best fit values in this case are 

H + AH = 67 km/s/Mpc, AH = 10 km/s/Mpc, r = 450 Mpc, Q = 0.29 . 

(3.15) 

Here the goodness of the fit is x 2 — 172.6 (x 2 /157 = 1.10). The confidence level 
contours with fi an d H fixed to their best fit values are displayed in Fig. |2l 




Figure 2: Confidence level contours in the LTB model with perfectly uniform present- 
day matter density: H {r) = 57 km/s/Mpc + AHe~ r/r °, tt M (r) = 0.29 (67 
km/s/Mpc) 2 /#o( r )- From OH . 



All these models have an inhomogenous Big Bang. One could also have an 
inhomogenous expansion with a spatially constant age of the universe by choosing 
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e.g. 



Hn(r) 



n M (r) 



H 



y/1 - Q M (r) - M (r)arsinh, I ^ImM 



;i - n M (r))y* 



^0 



1 + 5e~ r / r °) 2 



(3.16) 



The constraint of a simultaneous Big Bang leaves us with only one free function. The 
best fit values are OH 



H = 76.5 km/s/Mpc, 6 = 1.21, r = 1000 Mpc, Q c 



0.29 



(3.17) 



with x 2 = 175.5 (x 2 / 157 = 1 12 )- E q- <l3T7l) implies that the Hubble function 
if ( r ) varies from the value H (0) = 65 km/s/Mpc near us to its asymptotic value 
Hq{t r ) = 52 km/s/Mpc. The age of the universe is then t age = 1/H = 12.8 
Gyr. Similar values have also been found in the model of ref. [|28l . 

Thus simple and at least seemingly semirealistic LTB dust models can fit the 
supernova data. The point to note is that although the LTB equations of motion do not 
in general permit locally accelerated expansion, this does not exclude the possibility 
that there can be an effective, volume averaged acceleration, where a scale factor 
defined via the physical volume of some comoving region has a positive double time 
derivative [|43l . However, it can be shown there is no effective average acceleration 
[311 for the models considered above 6 . 



4 Towards more realistic LTB models 

Whether the supernova data combined with the CMB and LSS data would nevertheless 
require an accelerating universe is an open question; cosmological perturbation theory 
in LTB background is still non-existent. 

Some issues can be addressed, though. In particular, when the LTB metric models 
a local underdensity, one may assume that the evolution of perturbations is identical 
to that in a homogeneous universe until the time of last scattering. Adopting this 
approach, Alnes, Amarzguioui and Gr0n ||28T| have considered in an approximation 
constraints arising from the position 1% of first acoustic peak. They find a shift relative 
to the concordance ACDM model that is given by 

5 = ^ = 0.01419(1-00^, (4.1) 

where <1a is the angular diameter distance to the last scattering surface, given by Eq. 
(13.101) : this is the part that depends on the local underdensity, whereas r s , the sound 
horizon at recombination, and the (small) value of the parameter X can be obtained 



6 Although fitting the supernova data does not require accelerating expansion, for some profiles 
the LTB model may give rise to a suitably defined average acceleration B4l . For a discussion on 
backreaction in LTB models, see also l45ll . 
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Table 1: The best fit parameters of the locally underdense inhomogeneous model of 

El. 

from the conventional homogeneous model. To be in agreement with the WMAP 
observations, the shift parameter should be within the range S = 1.00 ± 0.01. The 
locally underdense model depends on the density contrast parameter Aa, functionally 
related to A(r, t) and specifying the difference between the two region, the transition 
point r from LTB to FRW, and the transition width Ar/r . A set of parameter values 
that yields a good fit both to the supernova data and the first CMB peak position can 
be found, as can be seen from Table CD Generically, for the void picture to work, one 
should have a local underdense region that extends at least up to the nearby supernovae 
or about 300-400 Mpc/h. 

These considerations hold if we occupy the exact center of the local LTB universe. 
For an observer that is located off-center, the universe appears to be anisotropic. 
Estimating the luminosity distance for an off-center observer is somewhat more 
complicated task than in the case of an observer at the center ll33l 146*1 . One finds an 
anisotropic relation between the redshifts and the luminosity distances of supernovae, 
which however yields only a mild constraint as up to about 20 % displacement from the 
center is consistent with the data [47 J. In contrast, the constraint from the CMB dipole 
appears to be very stringent, allowing only a displacement of about 15 Mpc from the 
center of the underdense bubble [33J. This result is obtained by assuming that all of 
the observed dipole a w ~ 10~ 3 is due to the displacement. A cancelation of the dipole 
due to our local peculiar motion towards the center of the underdensity is a possibility 
that would allow for a larger displacement. Whether such a peculiar motion can arise 
naturally or only by an accident, remains to be seen. 

For an off-center observer the direction towards the center of the bubble singles 
out a special axis. Therefore one could hope that a local LTB bubble could provide 
an explanation for the observed peculiar alignments of the CMB quadrupoles and 
octopoles Il48ll . Because of the smallness of the displacement allowed by the dipole, 
the quadru- and octopoles appear not to have enough power to explain their observed 
alignment [|33l . although again the conclusion depends on the assumption that our local 
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average motion has been accounted for correctly. 

Instead of a single underdensity, one could also consider an "onion model" with 
a homogeneous background density, on top of which there are density fluctuations 
which are periodic as a function of the radial coordinate. The observer sits in some 
generic position and looks at sources along the radial direction, and the LTB dust 
solution incorporates the entire Universe. To study this set-up, Biswas, Mansouri and 
Notari 11461 have derived an expression for the luminosity distance in an LTB metric 
for an off-centre observer. The corrections due to underdensities to light propagation 
were found to have a tendency to cancel far away from the observer because a radial 
light ray unavoidably meets both underdense and overdense structures. However, in 
the real universe light encounters hardly any structure, so the cancellations might 
be an artifact of the onion model. Since in the real universe the photon is mostly 
traversing voids it should get redshifted faster as the nonlinearities increase with 
time and thereby effectively produce an apparent acceleration. In the onion model 
one can nevertheless mimic an accelerating ACDM cosmology under certain special 
conditions: the observer has to be located around a minimum of the density contrast 
that is required to be quite high H6l . 

Yet another approach is the "Swiss cheese" model of the inhomogeneous universe, 
where each spherical void is described by the LTB metric. At the boundary of 
these regions the LTB metric is matched with the FRW metric that describes the 
evolution between the inhomogeneities. One can then seek for the modifications 
of the luminosity distance as the light passes through the underdense regions [|49l . 
In the extreme case where one assumes that light traverses the centers of all the 
inhomogeneities along its path, assuming that the locations of the source and the 
observer are random and inhomogeneities have sizes of order 10 Mpc, the relative 
increase of the luminosity distance is however just of the order of a few percent near 
z ~ 1. A qualitatively similar conclusion has been reached in [|50l . 

Structure formation and the smallness of CMB perturbations may in general pose a 
difficulty for LTB models. For instance, for a class of inhomogeneities a homogeneous 
universe is actually a late time attractor solution. This means that at earlier times matter 
density and/or Hubble rate tends to be even more inhomogeneous than today. Whether 
this presents an unsurmountable problem remains to be seen. Nevertheless, the LTB 
model serves as a reminder that the interpretation of the cosmological data is not only 
quantitatively but even qualitatively very much model dependendent. Therefore, all 
options should be carefully examined before firm conclusions can be drawn. This 
is true in particular for dark energy, which is both an observational and theoretical 
enigma. 
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